NAS A/CR- 1 999-208997 
ICASE Report No. 99-8 



The Dissipation Range in Rotating Turbulence 


Robert Rubinstein 
ICASE, Hampton, Virginia 

Ye Zhou 

ICASE, Hampton, Virginia 
and 

Tuskegee University, Tuskegee, Alabama 


Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, VA 

Operated by Universities Space Research Association 



National Aeronautics and 
Space Administration 

Langley Research Center 
Hampton, Virginia 23681-2199 


Prepared for Langley Research Center 
under Contract NAS 1-97046 


February 1999 


Available from the following: 


NASA Center for AeroSpace Information (CASI) 
7121 Standard Drive 
Hanover, MD 21076-1320 
(301)621-0390 


National Technical Information Service (NTIS) 
5285 Port Royal Road 
Springfield, VA 22161-2171 
(703) 487-4650 



THE DISSIPATION RANGE IN ROTATING TURBULENCE* 

ROBERT RUBINSTEIN^ AND YE ZHOU» 

Abstract. The dissipation range energy balance of the direct interaction approximation is applied to 
rotating turbulence when rotation effects persist well into the dissipation range. Assuming that RoRe l/ 2 << 

1 and that three- wave interactions are dominant, the dissipation range is found to be concentrated in the 
wavevector plane perpendicular to the rotation axis. This conclusion is consistent with previous analyses 
of inertial range energy transfer in rotating turbulence, which predict the accumulation of energy in those 
scales. 
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1. Introduction. Kraichnan [1] demonstrated that despite the analytical complexity of the equations 
of the direct interaction approximation, solutions representing the universal inertial and dissipation ranges of 
turbulence are readily accessible. In the inertial range, where energy transfer vanishes, the closure equations 
are satisfied by a power law spectrum with a constant energy flux from large to small scales of motion. 

In the dissipation range, viscous dissipation balances the energy input to each scale from nonlinear 
interaction. By assuming that the turbulent time scale is the viscous time scale 0 ~ (vk 2 )' 1 , the dissipation- 
range spectrum 

(1.1) E(k) ~ k 3 exp(-/3k/k d ) 

is obtained. In Eq. (1.1), k d = (e/i' 3 ) 1 '' 4 is the Kolmogorov scale; 0 will be used throughout to denote 
a universal constant, but not the same constant each time it appears. This spectrum is consistent with 
experimental and numerical data [2,3]. The theory also predicts that nonlinear interactions in the dissipation 
range are predominantly among nearly collinear wavevector triads; this conclusion is evaluated in [4, 5, 6]. 

Kraichnan [1] also proposed a near-dissipation range balance in which a nonlinear time scale replaces the 
viscous time scale. If this time scale is determined, consistently with Eulerian DIA by the sweeping hypothesis 
© ~ (V 0 k)~\ the energy spectrum E{k) ~ k 2 cxp{-k/k d ) is obtained. If the Kolmogorov time-scale is used 

instead, the result would be 

(1.2) £(fc) ~ /c 5/3 exp(-/c//c d ) 

Sirovich et al [7] propose an energy spectrum covering both the inertial and dissipation ranges which reduces 
to Eq. (1.2) in the dissipation range. 

2. Analysis. In rotating turbulence, the additional time-scale introduced by the Coriolis force makes a 
variety of dissipation range balances possible. Thus, although the very smallest scales will always be subject 
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to the viscous time-scale, leading to the spectrum of Eq. (1.1), a different result might be anticipated when 
wave interactions dominate nonlinearity well into the dissipation range. 

The standard elementary scaling arguments suggest when this condition applies. The wave frequency in 
rotating turbulence is determined by the dispersion relation of inertial waves as w(k) = Qk z /k. To begin, 
we will ignore the angular dependence of this frequency. Assuming the spectrum [8] for rapid rotation 
E(k) ~ (e:$2) 1 ' 2 k 2 and defining the rotational dissipation scale kd,u by the condition 

(2- 1 ) f d n vk 2 E{k)dk=e 

Jo 

there results k d> n ~ (eftlv 2 ) 1 / 2 . The viscous frequency scale exceeds the wave frequency scale at k' d Q defined 
b y n) 2 ~ Assume that this scale is much larger than k dy n so that 

/ x <-1/2 

( 2 * 2 ) k d t n/k r dti 1 ~ Q yi/2 ~ RoRe 1/2 << 1 

and consider the scales satisfying k d ,n < k < k dn . 

The direct interaction approximation (DIA) predicts the dissipation range energy balance [1] 

If r°° 

2vk Q(k) = - j dpdq <f(k - p - q) I dr x 

(2*8) Pimn (k)Q mr (p, T )Qns (<T t)G^' (k, T ) 

Standard notation is used in Eq. (2.3): ( k, r) = < u,(k, t)u 3 ( k', s) > /J(k+k') is the two-time correlation 

function, where r denotes the time difference r = t - s. The single-time correlation function is denoted by 
Qij (k), and the DIA response function is Gy(k, r). The tensorially isotropic form Q tJ (kj = Q(k)D i; fk) is 
assumed, where P«(k) = 6 tJ - k,k 3 k 2 and finally, P lmn (k) = fc m P in (k) + Jfc n P im (k). 

Eq. (2.3) is simply the general DIA energy balance, with eddy damping ignored in comparison to viscous 
damping. The assumption Eq. (2.2) means that wave turbulence theory [9] applies in the dissipation range. 
Therefore, the DIA response function Gij is simply the linear response 

( 2 - 4 ) G b ( k > T ) = {cos(2fiA: 2 T)/fc)£h(k) + sin(2fi* z r)/ifc)£°.(k)}i/(T) 

where H is the unit step function and the tensors ? are defined in terms of the Craya-Hemng basis 

e (1) (k) = kxfi/ | k x« | 

( 2 - 5 ) e< 2 >(k) = kx(kxll)/ | kx(kxft) | 

by 


( 2 . 6 ) 

Note that f 


£?. = ( 1 ) ( 2 ) _ ( 1 ) ( 2 ) 


= e t 0, ej 1) + e[ 


(2) J2) 


e- 


Fiji k). At lowest order in the weak turbulence approximation, 


( 2 *^) Q*i(k, t) — Gi m (k, r)Q mj (k) + Gj m (k, — r)Q mi (k) 

For the correlation function Q(k), we provisionally adopt Kraichnan’s dissipation range hypothesis [1] 
(2*8) Q{k) = K(k/kd,n) a e~^ k / kd ' n ) 
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where the dimensional constant k will be determined by the calculation. Substituting Eqs. (2.4)~(2.8) in the 
dissipation range balance Eq. (2.3), the right side will contain time integrals of the form 

/ oo p q z 

drexpirnij ± ± j} 

where the lower limit of integration insures that the result is real. We can assume [10] 

( 2 . 10 ) 


«<“• -■‘» = S«T ± 7 ± ^ 


As Eq. (2.10) shows, the integral in Eq. (2.3) vanishes unless the triad k,p,q satisfies the resonance 
condition 


( 2 . 11 ) 


^±^±^=0 

k p q 


But as Kraichnan [1] notes, the exponential ansatz Eq. (2.8) forces the triad to be nearly collinear: the 
quantities exp(-^fc) on the left side of Eq. (2.3) and exp -0{\ p | + | q |) on the right can be comparable 
only if | k |«| p | + | q j, whence the triangle inequality implies near collinearity of the vectors k,p,q. In 

this case, 


( 2 . 12 ) 

where zfcp ± q = k, and p • r = 0. therefore 
(2.13) 


p = 4- r 

k 


p i + 

r 2 

q l~ + -o) 


In view of Eq. (2.8). the integral in Eq. (2.3) decays exponentially with r. 

It is immediately evident that substitution of Eq. (2.12) into the resonance condition Eq. (2.11) yields 
the contradiction ±1 - 1 - 1 0 unless k. is nearly zero. Thus, the resonance and collinearity conditions 

can only be satisfied by nearh horizontal wavcvcctor triads. The right side of Eq. (2.3) is therefore nearly 
zero for non-horizontal wavrvivtors k. contradicting the assumption that the left side is isotropic. 

But by assuming instead that the energy spectrum is planar, so that instead of Eq. (2.8), 

(2.14) Qhk) = K{klk dt a) a e~^ k/kdn) 6(k z ) 
we find that Eq. (2.3) can 1« sati-da-d. since 

(2.15) J dp</q rt( k p - q)<5(p 2 )<5(q 2 ) = J ^ dp 6(p z )5(k z ) 

If the dissipation range excitation is confined to the plane k z = 0, the wave time-scale does not apply since 
the planar modes arc unaffected by rotation. The dissipation range balance must be computed assuming 
the viscous time-scale instead. With this understanding, the integrations in Eq. (2.3) are easily performed 
using the variables p, r in place of p, and lead to the result 


(2.16) 


K.k 2a+1 kj / v ~ vk 2 (k / kd,n) c 
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where the common exponential factor has been cancelled from both sides and a common factor of S(k z ) is 

cancelled using Eq. (2.15). The powers of k balance if a = 1, consequently k = v 2 . Therefore, the correlation 
function is 

( 2 - 17 ) <?(k) = Cv 2 {k/k dii i)e~ 0(k/kdn) 6{k z ) 

The effect of rotation is indirect: rotation excludes certain classes of vectors from the dissipation range, but 
the dissipation range dynamics in the remaining scales is independent of rotation. 

This conclusion was based on the impossibility of satisfying the three-wave resonance conditions by 
nearly collinear vectors, except when those vectors are nearly horizontal. In problems in which the dispersion 
relations do not permit three-wave interactions, perturbation theory leads to a modified equation of motion 
for waves with a higher order nonlinearity than the governing equations [11], If the governing equations 
are quadratically nonlinear, the result is a theory with a cubic nonlinearity. For such a theory of four-wave 
interactions, nontrivial resonances are always possible. 

It is therefore natural to ask whether four-wave interactions could play a role in the dissipation range of 
rotating turbulence. The possible role of four-wave interactions in the inertial range of rotating turbulence 

has been suggested by Yakhot [12], The dissipation range balance for four-wave interactions, which replaces 
Eq. (2.3) has the form 

o f f°° r°° r°° 

2iyk 2 Q tJ (k) = j dpdp'dqS(k - p - p' - q) / dt dt f ds x 

Pimn (k)P m ’rs(p)Pj'kl(k)Pk>p q (p)G 

mm' (p, t - s)G kk , (p, t' - s') x 

( 2 ‘ 18 ) G jr (k,t- t')Q rp (p',s - s')<? S9 (p - p',s - t)Q n i(q,t - s) 

in which only one representative term is written on the right side. The assumption of near collinearity takes 
the form 


p = ±£k + r 

k 

p ' = ±£k+ r ' 

( 2 - 19 ) q = ±fk-r 

k 

Quartets of this form can always be found to satisfy the four-wave resonance condition 

( 2 - 20 ) — ± — ±^± — = 0 

k p p' q 

Substituting the exponential ansatz Eq. (2.8) in the four-wave dissipation range balance Eq. (2.18), 

( 2 - 21 ) k 6 k 4 d , Q (k/k d , n ) 3a Sl- 2 K 2 = uk 2 (k/k d , n ) a 

Balancing the powers of k gives a = -2. The definition of k dM implies 

( 2 ‘ 22 ) Q{k) = C' D - 7 y /2 (k/k d , n )- 2 e- 0 ' k / k ^ 

In problems in which three- and four-wave interactions are both present, it can be assumed that the 
three- wave processes are dominant [13]. In this case, the conclusion that the dissipation range spectrum 
contains only horizontal vectors is consistent with the picture of energy transfer in rotating turbulence 
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proposed on the basis of closure studies by Cambon et al. [14, 15], namely that in the inertial range, energy 
is transferred to the wavenumber plane perpendicular to the rotation axis. Further theoretical support for 
this conclusion was provided by the instability principle of Waleffe [9]. 

The impossibility of dissipating energy in nearly vertical directions by energy transfer to smaller scales 
would naturally force the transfer of energy toward the horizontal plane, where viscous dissipation is possible. 
Although this picture of energy transfer in the inertial and dissipation ranges is self-consistent, it must be 
stressed that the conclusions about inertial range transfer do not require the strong restriction on the Rossby 
number required here. 

It should be noted that the prediction that the dissipation range is restricted to exactly horizontal vectors 
is a consequence of ignoring the angular dependence of the wave time-scale. In fact, since the wave frequency 
is Slk z /k, the condition k d ,a > k d occurs for wavevectors satisfying k z /k < cos 6 where cos 6 = er 1/2 /z^ 1/2 Q. 
For these modes, the viscous time-scale is dominant, and the three-wave dissipation range extends to the 
region -k cos 6 <k z <kcos8 instead of to the plane k z = 0. Note also that this argument shows that the 
scale in Eq. (2.17) should be replaced by the Kolmogorov scale k d . 

The approximate two-dimensionalization of the small scales of rotating turbulence under the limit defined 
by Eq. (2.2) raises the question of the relationship between the present results and the Taylor-Proudman the- 
orem [16], As noted by Smith and Waleffe [17], the large scales of rapidly rotating turbulence are always sub- 
ject to the Taylor-Proudman theorem since they are nearly steady. Whereas the Taylor-Proudman theorem 
requires the applicability of steady, linear dynamics, the present argument based on wave interactions allows 
both unsteadiness and nonlinearity. The possibility that a combination of large-scale two-dimensionalization 
due to the Taylor-Proudman theorem with two-dimensionalization of the small scales due to the impossibility 
of certain three- wave interactions leads to two-dimensionalization of all scales of motion in the extreme limit 
RoRe 1 / 2 << 1 is an interesting theoretical possibility which warrants further investigation. 
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